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Abstract 

This paper investigates the integration of on-line and off- 
line diagnostic algorithms for aircraft gas turbine engines. The 
on-line diagnostic algorithm is designed for in-flight fault 
detection. It continuously monitors engine outputs for 
anomalous signatures induced by faults. The off-line 
diagnostic algorithm is designed to track engine health 
degradation over the lifetime of an engine. It estimates engine 
health degradation periodically over the course of the engine’s 
life. The estimate generated by the off-line algorithm is used 
to “update” the on-line algorithm. Through this integration, 
the on-line algorithm becomes aware of engine health 
degradation, and its effectiveness to detect faults can be 
maintained while the engine continues to degrade. The benefit 
of this integration is investigated in a simulation environment 
using a nonlinear engine model. 

Introduction 

Early detection of faults is an important aspect in the health 
management of aircraft gas turbine engines. Undetected faults 
can lead the aircraft engine into an undesirable operating 
condition where compressor stall margin is reduced or turbine 
temperature is higher than the expected value. Under such a 
condition, safety and efficiency of engine operation may be 
compromised. Thus, it is critical to detect faults as early as 
possible and take the necessary corrective actions to avoid 
undesirable engine operation. 

To achieve a real-time fault detection capability, on-line 
diagnostic algorithms have been developed by several 
researchers (refs. 1 to 4). An on-line diagnostic algorithm is 
designed to run on an on-board engine computer in real time. 
It processes measured engine outputs to detect and, if possible, 
isolate a fault. Fault detection is pursued based on the fact that 
the measured engine outputs deviate from their reference 
condition values when an aircraft engine experiences a fault. 
An on-line diagnostic algorithm, however, encounters a 
challenge in achieving reliable performance. This challenge 
arises from the fact that the measured engine outputs are 
influenced not only by faults but also by engine health 
degradation. Engine health degradation is a normal aging 


process that all aircraft engines will experience due to usage 
and therefore is not considered as a fault, whereas a fault is an 
abnormal, unexpected event. Given the measured engine 
outputs, it is difficult to discern whether the engine output 
deviations are due to a fault or health degradation. Unless this 
issue is addressed, the on-line diagnostic algorithm will lose 
its effectiveness to detect faults as the engine degrades over its 
lifetime. 

One approach to address the above issue is to integrate the 
on-line diagnostic algorithm with an off-line trend monitoring 
algorithm. The task of an off-line trend monitoring algorithm 
is to track the engine health condition over the lifetime of an 
engine. It estimates the engine health condition based on 
steady-state engine output data recorded during flight (refs. 5 
to 7). Since health degradation is a gradual process, the off- 
line algorithm needs to update its estimate at a relatively low 
frequency, such as once per a few flights or days. This 
periodically updated knowledge of engine health condition can 
be used to adjust the on-line diagnostic algorithm. Through 
this integration, the on-line diagnostic algorithm becomes 
aware of health degradation, and its effectiveness to detect 
faults can be maintained while the engine continues to 
degrade. 

In this paper, an on-line fault detection algorithm and an 
off-line trend monitoring algorithm are integrated, and the 
benefit of this integration is investigated in a simulation 
environment. The on-line fault detection algorithm is based on 
the hybrid Kalman filter (refs. 8 and 9), and the off-line trend 
monitoring algorithm is based on the extended Kalman filter 
(ref. 10). In the following sections, the details of each 
algorithm are described, followed by the application of the 
design methodology to a large commercial aircraft engine 
model. Then, the effectiveness of the integrated diagnostic 
approach is evaluated using simulated examples of health 
degradation and faults. 

Nomenclature 

BST Booster 

EKF Extended Kalman filter 

F1KF Flybrid Kalman filter 

F1PC FTigh Pressure Compressor 
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HPT 

LPT 

OBEM 

P2 

P25 

Parab 

PLA 

PS3 

T2 

T25 

T3 

T49 

Tamb 

TCM 

TMHS23 

TMHS3 

TMHS41 

TMHS42 

TMHS5 

TMSHBC 

TMHSBL 

VBV 

VSV 

WF36 

WSSR 

XN12 

XN25 

XNH 

XNL 

e 

h 

Mcmd 

V 

X 

y 


High Pressure Turbine 
Low Pressure Turbine 
On-Board Engine Model 
Engine inlet pressure 
HPC inlet pressure 
Ambient pressure 
Power Lever Angle 
Combustor inlet static pressure 
Engine inlet temperature 
HPC inlet temperature 
Combustor inlet temperature 
LPT inlet temperature 
Ambient temperature 
Turbine clearance model 
BST metal temperature 
HPC metal temperature 
HPT nozzle metal temperature 
HPT metal temperature 
LPT metal temperature 
Combustor case metal temperature 
Combustor liner metal temperature 
Variable bleed valve 
Variable stator vane 
Fuel flow 

Weighted Sum of Squared Residuals 

Fan speed, measured 

Core speed, measured 

Core speed, actual 

Fan speed, actual 

Environmental parameter vector 

Health parameter vector 

Control command vector 

Sensor noise vector 

State variable vector 

Sensor output vector (controls/diagnostics) 
Sensor output vector (ambient/engine inlet) 


Problem Background: Influence of Faults 
And Health Degradation 

Although aircraft gas turbine engines are highly reliable, 
they may encounter numerous types of faults during flight. 
Sensor faults are most likely to occur (ref. 1), but rotating 
components and actuators may also experience faults. When 
an aircraft engine experiences a fault, the measured engine 
outputs deviate from their reference condition values. 
Therefore, fault detection can be accomplished by monitoring 
for shifts in measured engine outputs induced by faults. This 
is, however, a challenging problem since shifts in measured 
engine outputs are induced not only by faults but also by 
engine health degradation. Engine health degradation is the 
result of usage and manifests itself as degraded performance 
of engine components such as compressors and turbines. As it 
worsens gradually over time, the deviation of the measured 
engine outputs from their reference condition values increases. 


Figure 1 shows the influence of health degradation and 
faults on a measured engine output. The horizontal axis 
indicates the sample index with a relatively long time interval 
(e.g., a few flights or a few days). The vertical axis indicates 
the measured engine output at a specific flight condition. The 
reference is an expected value of that engine output, and it is 
derived at a specific health baseline (e.g., nominal health 
condition). The gradual deviation of the engine output from 
the reference indicates the influence of health degradation. 
The discrete jump from point '‘a” to point “b” at the k th sample 
indicates fault occurrence. At the sample point k, the engine 
output shift Aic from the reference is induced by health 
degradation while a fault induces additional shift 5A. Thus, the 
total shift (A k + 5A) is observed through a sensor. An on-line 
diagnostic algorithm, in general, has to process this total shift 
in order to detect the fault. If the algorithm does not have the 
capability to handle the health degradation-induced shifts in 
the sensor measurements, it may either miss the fault or 
generate a false alarm. Missed detection may occur when the 
fault is masked by health degradation. A false alarm may be 
generated when health degradation is misdiagnosed as a fault. 
Either way, the on-line algorithm will eventually lose its 
diagnostic effectiveness. 

To address the problem described above, the integration of 
an on-line fault detection algorithm and an off-line trend 
monitoring algorithm was proposed in references 8 and 9. The 
integration architecture of the two diagnostic algorithms is 
shown in figure 2. The on-line algorithm processes measured 
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Figure 1. — Influence of health degradation and faults. 
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Figure 2. — Architecture of integrated diagnostic approach. 
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engine outputs in real time to detect faults during flight. Such 
an algorithm is generally developed at a specific health 
baseline, and its reference for measured engine outputs is 
fixed. Therefore, the measurements processed in real time 
contain health degradation-induced shifts, which gradually 
increase over the lifetime of an engine. To account for the 
influence of health degradation, the health baseline of the on- 
line algorithm must be updated periodically so that the 
algorithm can operate in the vicinity of the degraded engine’s 
health condition. Through the health baseline update, the 
reference for measured engine outputs is shifted to the vicinity 
of the degraded engine’s outputs (for instance, the reference is 
shifted to the vicinity of point “a” in fig. 1). Consequently, the 
health degradation induced-shifts that the on-line algorithm 
encounters are significantly reduced as compared to the fixed 
reference case. 

The health baseline which is needed to update the on-line 
algorithm is provided by the off-line trend monitoring 
algorithm as shown in figure 2. The off-line algorithm 
processes data post-flight to estimate the health condition of 
the degraded engine. This data may contain engine outputs 
recorded at a single or multiple operating conditions. Since 
health degradation is a gradual process, the frequency at which 
the off-line algorithm operates is much lower than that of the 
on-line algorithm. 

The algorithms used in this paper for on-line and off-line 
diagnostics are based on model-based estimation techniques. 
The on-line algorithm for fault detection is based on the 
hybrid Kalman filter (HKF) technique. The HKF is a uniquely 
structured Kalman filter which lends itself to the health 
baseline update. The off-line algorithm for trend monitoring is 
based on a derivative of the extended Kalman filter (EKF) 
technique. A description of each algorithm is given in the 
following sections. 

Integration of On-Line and Off-Line 
Diagnostic Algorithms 

The on-line and off-line diagnostic algorithms discussed in 
this section are based on a well-known estimation technique: 
the Kalman filter. Both algorithms update their estimations as 
new measurement data become available. The major factor 
that differentiates these algorithms is the frequency at which 
the measurement data is sampled. Because of this difference, 
the aircraft engine is described in a different form for each 
algorithm. In this section, the description of the aircraft engine 
is given first. Then, the descriptions of on-line and off-line 
diagnostic algorithms are given. 

Aircraft Engine Description 

An aircraft engine under consideration for on-line 
diagnostics is described by nonlinear equations of the 
following form: 


X ~ f{x,h k , W and , ^) 

y = g(x, h k ,u cmd ,e) + v 

The vectors x, u cmd , and e contain state variables, control 
command inputs, and environmental parameters, respectively. 
Health degradation of the engine is represented by the vector h 
which contains health parameters. Health parameters are 
efficiencies and flow capacities of the engine components 
such as compressors and turbines. As they deviate from the 
nominal (healthy) baseline, the performance delivered by each 
component degrades. Since health degradation progresses 
gradually over time, the health parameter vector is treated as a 
constant at the time scale where on-line diagnostics is 
performed. The subscript k of the health parameter vector 
indicates the sample index at the time scale where off-line 
diagnostics is performed. For given inputs, the nonlinear 
functions / and g generate the state derivative vector x and 
the sensor output vector y, respectively. The sensor outputs are 
corrupted by the white noise vector v. 

While the aircraft engine is described in the form of 
equation (1) for on-line diagnostics, it is described in the 
following form for off-line diagnostics: 

Y k =gss{hk,U k ,E k )+V k (2) 

The vectors Y, U, and E contain steady-state values of sensor 
outputs, control command inputs, and environmental 
parameters, respectively. The subscript k indicates the sample 
index. The nonlinear function g ss is equivalent to g in 
equation (1) at steady state. The vector V represents “steady- 
state” white noise, and its magnitude is much smaller than that 
of actual noise v contained in the raw measurements. The 
steady-state noise represents small biases that exist at the time 
of data sampling. 

In the off-line trend monitoring problem, the health 
parameter vector is treated as time variant; the health 
condition can change over the time interval between the 
sample points. Because of the gradual nature of health 
degradation, the time interval can be rather long. Depending 
on how the off-line trend monitoring is performed, the time 
interval can be a few flights or a few days. 

On-Line Algorithm: Hybrid Kalman Filter 

The objective of on-line fault detection is to detect faults as 
early as possible from the observed engine outputs. To achieve 
this objective, the sensor outputs in equation (1) are 
continuously processed in real time so that any shift induced 
by faults can be detected. The on-line diagnostic algorithm, 
which pursues the above objective, is based on the hybrid 
Kalman fdter (HKF) technique (refs. 8 and 9). The HKF is a 
uniquely structured Kalman filter which utilizes a nonlinear 
on-board engine model (OBEM) of the following form: 
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*OBEM = /(*OBEM , ^OBEM , u cmd , z ) (3) 

TOBEM = g'(^OBEM > ^OBEM > Ucmd > z ) 

This form is similar to the assumed form of the aircraft engine 
in equation (1). The vector xqbem contains the state variables 
of the OBEM while the vector z contains the measured 
parameters which define the flight condition. The vector 
/7 0 bem represents the health baseline of the OBEM. By 
utilizing the state and output vectors of the OBEM, the 
following E1KF equation is obtained: 

x = A(x-x OBEM )+K{y-y) ^ 

y = CyX — -T 0 BEM ) + T OBEM 

The vectors x and y represent the estimate of x and y of the 
plant in equation (1), respectively. The state-space matrices A 
and C are derived through the linearization of the nonlinear 
plant model (in the present case, OBEM). The matrix K 
represents the Kalman gain. Similar to the general linear 
Kalman filter approach, the matrices A, C, and K are derived 
at multiple operating conditions to cover a wide operating 
range of the aircraft engine. Once derived, these matrices 
are saved in table lookup form for real-time execution of 
equation (4). 

The Kalman filter, in general, is able to accurately estimate 
the sensor outputs as long as the plant operates under the 
conditions expected by the Kalman filter. When the plant 
experiences unexpected events such as faults, the estimation 
performance of the Kalman filter becomes poor, causing an 
increase in the residuals, [y - y) . Therefore, fault detection 
can be achieved by detecting such an increase in the residuals. 
To detect fault-induced shifts in the residuals, weighted sum 
of squared residuals (WSSR) signals are constructed as 
follows: 

WSSRhkf ={y-y) T Z- l {y-y) 

WSSR obem ={y- Tobem ) r ^ -1 (y - jobem ) 

where 

X = diag[o 2 ] 

The vector a represents the standard deviation of the sensor 
measurements. The square matrix X normalizes the residual 
vectors (y - y) and (t-Jobem) • The fault indicator signals, 
WSSRhkf and WSSRobem, indicate the existence of faults; 
they will increase when the engine experiences faults. To 
detect faults, these signals are compared against pre- 
established thresholds. Fault detection is declared when at 
least one of the two fault indicator signals exceeds a threshold. 
As demonstrated in reference 9, the utilization of the two fault 


indicator signals improves the fault detection capability since 
these signals have different sensitivities to different types of 
faults. 

The EIKF designed at any particular health condition 
eventually loses its estimation accuracy as the engine degrades 
over its lifetime, and thus the fault indicator signals will 
increase. This is due to the fact that the discrepancy between 
the design health condition of the EIKF and the actual health 
condition of the engine increases over time. In the HKF 
approach, the design health condition is defined by the vector 
/?obem at which the OBEM operates and also by the health 
condition at which the matrices A, C, and K are derived. 

In references 8 and 9, it is shown that the HKF is able to 
maintain its estimation performance as long as the health 
baseline of the OBEM remains in the vicinity of the engine’s 
health condition. Thus, the HKF -based on-line fault detection 
algorithm is able to maintain its effectiveness by periodically 
updating the health baseline of the OBEM to a new condition 
as follows: 

hoBEM = hk (6) 

The vector h k represents the estimate of the health parameter 
vector h k in equation (1). This estimate is provided by the off- 
line trend monitoring algorithm, which will be discussed in the 
following section. The update process for the HKF is 
relatively simple; it is completed by feeding the estimated 
health degradation values into the OBEM. The matrices A, C, 
and K do not need to be redesigned as the health baseline of 
the OBEM is updated. This relatively simple update process is 
a major benefit of the HKF approach. 

Off-Line Algorithm: Extended Kalman Filter 

The objective of off-line trend monitoring is to track the 
health condition of an engine over its lifetime. This objective 
is accomplished by accurately estimating the health parameter 
vector hk in equation (2) from the steady-state sensor outputs 
Y k . The off-line diagnostic algorithm used in this study is 
based on a derivative of the extended Kalman filter (EKF) 
technique (ref. 10). The EKF is similar to the general linear 
Kalman filter except that a nonlinear plant model is 
successively re-linearized so that the validity of a linear plant 
representation is maintained as the plant moves from one 
operating condition to another. Therefore, a nonlinear plant 
model functions as part of the algorithm in the EKF approach. 
The nonlinear aircraft engine model used in the current study 
is given in the following form: 

Yk=gss[h,U k ,Z k ) (7) 

This form is similar to the assumed form of the aircraft engine 
in equation (2). The vector U is the same vector that appears in 
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equation (2). The vector Z contains the steady-state values of 
measured parameters that indicate the flight condition. The 

subscript k indicates the sample index. The vector h 
represents the estimate of the health parameter vector h, and 
this vector must be derived by the off-line trend monitoring 
algorithm. 

The above nonlinear plant model can be presented in the 
following alternative form: 

Yk = G k \ \k-\ [h - h- 1 )+ Yk\k-1 ( 8 ) 

where 

Yk\k-\ =8ss ( /? A--t ,U k ,Z k ) (9) 

G k\k-1 = — gss[h-i,Uk,Z k ) (10) 

dh 

The vector h k _ i represents the health parameter vector 
estimated at the previous sample point, k— 1. The matrix 
G k\k-\ represents the influence coefficient matrix derived 

from the nonlinear function, g ss (h k _ l ,U k ,Z k )• Given the 

vectors h k _ \ , Uk , Z k , and Y k , the health parameter vector at 
the sampling point k is estimated as follows (ref. 6): 

h k =h k -i+L^k-Yk\k-x) (11) 

where 

<i2) 

The matrices Pq and R represent the error covariance and 
measurement covariance, respectively. In the general Kalman 
filter approach (both linear and extended), the error covariance 
matrix is updated as the new measurement data become 
available (ref. 10). This matrix, however, may diverge unless 
the number of available measurements is greater than or equal 
to the number of health parameters to be estimated. In general, 
the number of sensors installed on an aircraft engine is less 
than the number of health parameters. To avoid the 
divergence, the error covariance matrix is kept constant. 
Otheiwise, the matrix update must be done with a special 
treatment, such as the soft-constraint approach demonstrated 
in reference 1 1 . 

Equations (9) through (12) represent the EKF algorithm. 
The algorithm updates its estimate as new data Y k , Uk, and Z k 
become available. To check the validity of the estimated 
health condition, the following weighted sum of squared 
residuals is computed: 


WSSR ekf = (y k - Y k f 2 -1 (y k -Y k ) ( 13 ) 

This signal is compared against a pre-determined threshold. If 
the signal exceeds the threshold, it is considered that the 
estimated health condition is inaccurate, and thus further 
investigation is needed to determine a cause of inaccuracy 

After validation, the estimated health condition h k is used 
to update the health baseline of the OBEM as described by 
equation (6). Through this update, the effectiveness of the on- 
line fault detection algorithm is maintained in the presence of 
health degradation. The estimated health condition also gives 
valuable information for engine health management. The 
snapshot estimation at a sample point can be used to consider 
whether a maintenance action is necessary or not. Thus, 
condition-based maintenance is possible instead of time -based 
maintenance. Engine health management can be further 
improved by interpreting the history of the estimated health 
condition. The estimated health condition from the past to the 
present reveals the degradation profile under which a specific 
engine operated. From this information, the health condition at 
a future point can be anticipated, and necessary actions may be 
taken to improve the safety or efficiency of the engine 
operation. 

Application of the Integrated Methodology 
to an Aircraft Engine Model 

The algorithms described in the previous sections require a 
nonlinear plant model. More specifically, the on-line 
algorithm requires a nonlinear model that can be executed in 
real time. In this section, a description of the aircraft engine 
model is given. For the results presented in this paper, the 
same model used by the diagnostic algorithms is also used to 
represent a real engine. The health degradation profile that the 
real engine may undergo is also discussed. 

Engine Model 

The engine model used in this paper is a nonlinear 
simulation of an advanced high-bypass turbofan engine, a 
typical power plant for a large commercial aircraft. This 
engine model has been constructed as a component level 
model (CLM), which consists of the major components of an 
aircraft engine. The CLM represents highly complex engine 
physics while being designed to run in real time. Engine 
performance deviations from the nominal health baseline are 
modeled by adjustments to efficiency and flow capacity 
scalars of the following five components: Fan (FAN), Booster 
(BST), High-Pressure Compressor (HPC), High-Pressure 
Turbine (HPT), and Low-Pressure Turbine (LPT). There are a 
total of 10 adjustments that are called health parameters. The 
engine state variables, health parameters, actuator variables, 
and environmental parameters are shown in Table 1 . 
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TABLE 1.— ENGINE MODEL VARIABLES 


State Variables 

XNL, XNH, TMHS23, TMHS3 
TMHSBL, TMHSBC, TMHS41 
TMHS42, TMHS5 

Health Parameters 

FAN efficiency, FAN flow capacity 
BST efficiency, BST flow capacity 
HPC efficiency, HPC flow capacity 
HPT efficiency, HPT flow capacity 
LPT efficiency, LPT flow capacity 

Actuators 

WF36, VBV, VSV 

Environmental Parameters 

Altitude, Mach Number 
Ambient Temperature 


There are a total of 1 1 measured parameters (v and z) that 
are available to the digital control unit of this engine. Table 2 
shows seven sensors (y) along with their standard deviations 
given in percent of steady-state values at the ground maximum 
power condition. The control actions and diagnostics are based 
on those sensed variables. Table 3 shows four additional 
measured parameters (z) along with their standard deviations 
given in their actual engineering units. These four parameters 
indicate the ambient and engine inlet conditions. The 
measurements of the inlet condition, T2 and P2, are used for 
parameter corrections (ref. 12). 

TABLE 2.— STANDARD DEVIATIONS OF 
CONTROLS AND DIAGNOSTICS SENSORS 


(a in Percent of Steady-State Values at 
Ground Maximum Power Condition) 


Sensors (y) 

a (%) 

XN12 

0.25 

XN25 

0.25 

P25 

0.50 

T25 

0.75 

PS3 

0.50 

T3 

0.75 

T49 

0.75 


TABLE 3.— STANDARD DEVIATIONS OF 
AMBIENT AND ENGINE INLET SENSORS 


(o in Actual Units) 


Sensors ( z ) 

a 

Tamb 

5.0 °F 

p 

1 amb 

0.1 psi 

T2 

5.0 °F 

P2 

0.1 psi 


control system adjusts three actuation variables to cause the 
corrected measured fan speed to match the desired value. The 
closed-loop system runs in a simulation environment at the 
frequency of 50 Hz. 

The two engine models used for on-line and off-line 
diagnostics are identical except that they operate at different time 
scales. These two engine models operate at estimated health 
conditions, and their flight condition is specified by three 
measured parameters: T amb , P lim b, and T2. From these three 
measurements, the engine models calculate the altitude, Mach 
number, and the temperature deviation from the standard day 
condition. The two engine models receive the three control 
commands (Table 1) generated by the control system. 

In the current study, the same engine model is used to 
represent a real engine and the engine model for diagnostics. 
This, however, does not represent a realistic application 
environment; there will typically be a mismatch between the real 
engine and the engine model due to un-modeled and incorrectly 
modeled physical phenomena. To represent such a mismatch, a 
turbine clearance model (TCM) is added only to the engine 
model representing a real engine. The TCM captures the turbine 
clearance behavior with high fidelity, and its addition results in 
output mismatch between the real engine and the engine model 
at steady-state conditions and during transients. 


Health Degradation Profile 


Each aircraft engine degrades differently depending on the 
environmental conditions the engine is exposed to. For instance, 
a high concentration of sand or salt can accelerate fan or 
compressor degradation while tropical heat can accelerate 
turbine degradation. Thus, there is no such thing as a typical 
degradation profile that engines will undergo (refs. 13 and 14). 
To capture this randomness, health degradation profiles 
considered in this paper are generated as follows: 


K = h 


i 

k - 1 



if k> 0 


K= 8 


i 

k 9 


if k = Q 


(14) 


where 


k = 0, 1, • • • , q 
i =l,2,---,p 


The nonlinear engine model is used in the subsequent 
sections to represent a real engine in equation (1) and the 
engine models in equations (3) and (7) for on-line and off-line 
diagnostics. The engine model representing a real engine 
operates at given health conditions, and its flight condition is 
specified by the three environmental parameters listed in 
Table 1. This engine model operates in closed loop with a 
control system described in reference 9. In the current control 
architecture, the power lever angle (PLA) is converted to 
desired corrected fan speed (an indicator of thrust). The 


The integers q and p indicate the number of sample points and 
the number of health parameters, respectively. Equation (14) 
gives the value of the i th health parameter at the k! h sample point. 
The scalar if specifies the maximum amount that the i th health 
parameter can deviate over q sample points. The random number 
e is uniformly distributed over the range between 0 and 1 . At the 
initial point k= 0, each health parameter is specified by 5, a 
random number uniformly distributed over the range of ±0.5%. 
The initial health condition represents engine-to-engine 
variations due to manufacturing tolerance. 
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With the above equation, each health parameter undergoes a 
unique degradation profile. The degradation profile depends on 
the scalar q specified for each health parameter and the random 
number s. 

Performance Evaluation 

In this section, the perfonnance of the off-line trend 
monitoring algorithm is investigated first, using a specific 
degradation profile. The off-line algorithm estimates the engine 
health condition at each sample point. The estimated health 
condition is then used to update the on-line fault detection 
algorithm. After the health baseline update, the on-line 
algorithm’s capability to detect sensor faults is evaluated. The 
engine to be diagnosed by both algorithms contains the TCM in 
order to represent model-plant mismatch that exists in the real 
environment. 

Off-Line Trend Monitoring 


does not exist. The trend monitoring algorithm estimates the 
health parameters in such a way that the difference between the 
engine and model outputs is minimized in a least squares sense, 
while simultaneously minimizing the difference in the estimated 
health conditions at the current and previous sample points. 1 This 
approach may produce health parameter estimates that are 
inaccurate. To work around this problem, the following 
constraint is imposed: the estimated health parameters are 
constrained to change in the specific direction expected from the 
heuristic knowledge regarding health degradation. All estimated 
health parameters are not allowed to increase from the estimation 
values at the previous sample point, except for HPT flow and 
LPT flow which are not allowed to decrease. If an estimated 
health parameter changes in the direction opposite of that 
expected, then the estimation value at the previous sample point 
is used as the current estimation value. 

Figure 3 shows the actual and estimated health conditions 
over 300 sample points of engine life. The estimated health 
parameters were generated under the constraint previously 
discussed, except for the first five sample points. During the 


The off-line trend monitoring algorithm in equations (9) 
through (12) receives the steady-state vectors Y k , U k , and Z k and 
estimates the engine health condition at each sample point. The 
time interval between the sample points is rather long as 
discussed earlier. The steady-state data used in this section are 
generated from the engine model in equation (1) through the 
following steps. First, a degradation profile that the engine model 
will undergo over its lifetime ( q sample points) is generated 
using equation (14). The degradation profile specifies a unique 
health condition at each sample point. Then, at each sample 
point, the engine model is ran at a steady-state cruise condition. 
The operating condition of the engine is specified by the health 
condition, PLA, and environmental parameters in Table 1. From 
the steady-state simulation run, one-second time histories of y, 
u cm j, and z are collected. Then, the noise-corrupted time histories 
are averaged to obtain the steady-state vectors Y k , U k , and Z k . 

In the real environment, the steady-state flight condition of the 
engine can vary from flight to flight even at cruise. To capture 
this variation, PLA, altitude, Mach number, and ambient 
temperature are randomly varied from sample point to sample 
point. The PLA is varied over the range between 62° and 66° 
(intermediate power at cruise). Altitude is varied over the range 
between 33,000 and 37,000 ft. Mach number is varied over the 
range between Mach 0.83 and 0.85. The ambient temperature is 
varied within ±30 °F from the standard day condition. 

During a preliminary evaluation study, it was found that the 
off-line trend monitoring algorithm had difficulty in estimating 
FAN flow and LPT flow. When the actual FAN flow and LPT 
flow deviated toward a certain direction from their initial 
condition values, the estimated FAN flow and LPT flow did not 
follow the actual values. Instead, these estimated values deviated 
toward the direction opposite of the actual FAN flow and LPT 
flow. In the current estimation problem where the number of 
unknown parameters (health parameters) is greater than the 
number of known parameters (sensor outputs), a unique solution 
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Figure 3. — Actual and estimated health 
degradation over the engine’s lifetime. 


‘The algorithm given by equations (11) and (12) is equivalent to the 
maximum a posteriori estimator (refs. 15 and 16) and the minimum 
variance estimator (refs. 15 and 17). 
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first five sample points, the health condition of the engine is 
maintained at the initial health condition representing engine- 
to-engine variation due to manufacturing tolerance. The 
estimated health parameters which are initialized at the 
nominal health condition are allowed to change in any 
direction; constraints are not imposed on the estimated values 
during this period. The first five samples are simply used to let 
the estimated health parameters settle to any values in the 
presence of initial model-plant mismatch, which not only 
includes health condition mismatch but also the mismatch due 
to the TCM. 

As can be seen in figure 3, the algorithm is able to follow 
the degradation profile quite well for three components: BST, 
HPC, and HPT. There are estimation errors, but the errors 
remain almost constant across the sample points. Thus, the 
deviation from the initial condition can be accurately tracked 
for these three components. The algorithm displays some 
difficulty in estimating the health parameters of FAN and LPT 
components; the estimation errors increase with time. The 
cause of this divergence is not clear. It may be due to the 
underdetermined nature of the current problem, or the location 
of the sensors may not be suited for estimating these health 
parameters. When this study was repeated using different 
patterns of health degradation profile, a similar tendency was 
observed. A total of five patterns of degradation profile, 
generated from equation (14), was investigated. 

Health Baseline Update of On-Line Fault Detection 
Algorithm 

The on-line fault detection algorithm in equations (3) and 
(4) receives the vectors y, u cm </, and z and generates the sensor 
output estimates and the on-board model outputs at each 
sample point with the time increment of 0.02 sec (50 Hz). 
Based on the residual vectors (y - y) and (y - _y OBEM ) , the two 
fault indicator signals in equation (5) are computed, and they 
are compared against pre-established thresholds. If at least one 
of the two fault indicator signals exceeds the threshold for 25 
consecutive time steps, fault detection is declared as 
demonstrated in reference 9. 

Before evaluating the fault detection capability of the on- 
line algorithm, the benefit through the integration of the on- 
line and off-line algorithms is investigated. As discussed 
earlier, the on-line fault detection algorithm will eventually 
lose its effectiveness as the engine degrades over time. This 
happens because the fault indicator signal values increase with 
the increase in discrepancy between the engine and the 
OBEM. To maintain diagnostic effectiveness, the health 
baseline of the OBEM must be updated periodically. This 
update is completed by feeding the estimated health condition 
values from the off-line trend monitoring algorithm to the 
OBEM. The importance of this health baseline update is 
illustrated in figure 4 where, over its lifetime of 300 sample 
points, the aircraft engine is subject to the specific health 
degradation profile shown in figure 3. At each sample point, 
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Figure 4. — Influence of engine health degradation on steady- 
state values of WSSR h «f and WSSRqbem- 


the engine and the on-line diagnostic algorithm are run at a 
specific cruise condition (35,000 ft altitude, Mach 0.84, 65° 
PLA). Figure 4 shows the steady-state values of the fault 
indicator signals for the following two cases: 1) the health 
baseline of the OBEM maintained at the nominal condition 
and 2) the health baseline of the OBEM updated to the 
estimated health parameters shown in figure 3. The dashed 
lines indicate fault detection thresholds. The threshold values 
were determined using health condition mismatches between 
the engine and the OBEM as done in references 8 and 9. 

As can be seen in the figure, if the health baseline of the 
OBEM is not updated from the nominal condition, WSSRhkf 
exceeds the threshold around the 130 th sample point while 
WSSRobem exceeds the threshold around the 30 th sample 
point. This means that the fault detection algorithm will start 
misdiagnosing health degradation as a fault around the 30 th 
sample point. To avoid generating false alarms, the thresholds 
must be set at much higher values. However, doing so will 
compromise the fault detection capability and will increase the 
missed detection rate when the engine is at or near the nominal 
health condition. 

When the health baseline of the OBEM is updated 
periodically, the fault indicator signals remain fairly constant 
over the lifetime of the degrading engine. In figure 4, the 
OBEM is subject to the estimated health condition shown in 
figure 3, therefore, the estimation error in some of the health 
parameters increases over time. This increase in the estimation 
error seems to have a very minor influence on the fault 
indicator signals. The reason for this may be attributed to the 
nature of the off-line trend monitoring algorithm. The off-line 
trend monitoring algorithm estimates the health parameters in 
such a way that the difference between the engine and model 
outputs is minimized in a least squares sense (refs. 15 to 17). 
Thus, the estimated health parameters are not necessarily 
accurate, but they are ones that minimize the difference 
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between the engine and model outputs. When these estimated 
health parameters are fed into the OBEM, the outputs of the 
OBEM match well with the engine outputs. Consequently, the 
fault indicator signals remain small. 

Updating the health baseline of the OBEM has the effect of 
shifting the engine output reference to a new value. As 
illustrated in figure 1, the original reference at which the 
algorithm is designed becomes obsolete as the engine 
continues to degrade. By shifting the reference to a new value, 
such as point “a” in the figure, the on-line algorithm can detect 
the fault-induced shift, 5A, more reliably. 

On-Line Fault Detection: Sensor Bias Case 

The on-line algorithm’s capability to detect faults is 
evaluated in this section using sensor biases. For this 
evaluation, the closed-loop engine is trimmed at a specific 
cruise condition (35,000 ft altitude, Mach 0.84, 65° PLA) in 
the presence of health degradation and also a bias in a single 
sensor. Then, the engine and the on-line algorithm are ran for 
100 sec at steady state. The smallest magnitude bias that is 
detected for each sensor is shown in Table 4. In this study, the 
aircraft engine and the OBEM are subject to, respectively, the 
actual and estimated health conditions in figure 3 at three 
sample points: £=100, 200, and 300. 

From Table 4, it can be observed that the on-line algorithm 
performs consistently in the presence of health degradation; 
the detected bias magnitudes are almost the same at three 
health conditions, except for slight variations in the XN12 and 
XN25 sensors. For some of the sensors, detected bias 
magnitudes are not symmetric between the positive and 
negative directions. As discussed in reference 8, health 
condition mismatch between the engine and the OBEM has an 
effect of counteracting a sensor bias to some extent. As a 
result, depending on the bias direction, it takes smaller or 
larger magnitude for the fault indicator signals to exceed the 
fault detection thresholds. This lack of symmetry, however, is 
not as prominent as the result in reference 8. 


TABLE 4.— MINIMUM SENSOR BIAS 
DETECTED AT CRUISE WITH 65° PLA 



£ = 100 

£=200 

£ = 300 


(<fl 

(o) 

(a) 

XN12 

2.6 / -4.5 

2.9/ -4.5 

2.9/ -4.4 

XN25 

3.2/ -3.8 

3.4/ -3 .7 

3.3/ -3.6 

P25 

1.3/-1.3 

1.3/-1.3 

1.3 / -1.3 

T25 

2.1 1-22 

2.0/ -2.1 

2.0/ -2.1 

PS3 

1.7 / —1.6 

1.7/-1.7 

1.7 / — 1.7 

T3 

2.6/ -2.6 

2.6/ -2.6 

2.6/ -2.6 

T49 

2.9/ -3.5 

2.9/ -3.5 

2.9/ -3.5 





Tamb 

0.3/ -0.3 

0.3/ -0.3 

0.3/ -0.3 

P 

r amb 

0.5/ -0.5 

0.5/ -0.5 

0.5/ -0.5 

T2 

0.4/ -0.4 

0.4/ -0.4 

0.4/ -0.4 

P2 

/ — 

/ 

— / 


As mentioned earlier, three sensor measurements, T arnh , 
Pamb, and T2 are used to define the flight condition for the 
OBEM. A bias in these sensors, therefore, causes the OBEM 
to operate at a flight condition different from the true 
condition. The table shows that biases of small magnitude are 
detected for these three sensors. The on-line algorithm is not 
able to detect a bias in the P2 sensor. This sensor is currently 
used for the correction of pressure measurements. Since both 
measured and estimated pressure values are corrected by the 
P2 value, a bias in this sensor does not increase the residuals. 
The reader should be reminded that the fault detection 
algorithm will not identify the biased sensor; it only indicates 
that a fault exists. After the detection of a fault, the identity or 
severity of the detected fault must be classified through the 
fault isolation process (ref. 18). 

Discussion 

The study in this paper reveals the benefit of integrating the 
on-line and off-line diagnostic algorithms. Through this 
integration, the on-line fault detection algorithm is able to 
maintain its diagnostic effectiveness as the aircraft engine 
degrades over its lifetime. The integration of the two 
algorithms, however, is not a trivial step. In this section, the 
practical aspects of the integrated approach are discussed. 

One issue that influences the practicality of the integrated 
approach is the complexity involved in the update process of 
the on-line algorithm. If the update process involves a re- 
design of the entire on-line algorithm, such as the case for the 
piecewise linear Kalman filter approach (ref. 18), the re- 
design process may have to take place on a ground-based 
computer with human intervention. Furthermore, after the re- 
design process, the on-line algorithm must be uploaded to an 
on-board engine computer. Such complexity makes that 
approach impractical. It is desirable to keep the update process 
as simple as possible. In the case of the HKF, the algorithm is 
updated through a relatively simple process: by feeding the 
estimated health condition values to the OBEM. Thus, the 
HKF is well suited for the integrated approach. 

Another issue that influences the practicality of the 
integrated approach is the location of the off-line trend 
monitoring algorithm. In general, a trend monitoring algorithm 
resides on a ground-based computer. In that case, the data for 
post-flight analysis and the off-line diagnostic results must be 
transmitted between the ground-based and on-board 
computers. There would be a number of communication and 
certification issues that would have to be addressed to achieve 
this integration. 

The difficulties associated with off-board to on-board data 
transmission can be avoided if the off-line trend monitoring 
algorithm is capable of running on an on-board engine 
computer. Such a capability heavily depends on the off-line 
algorithm’s robustness. As discussed earlier, trend monitoring 
is generally an underdetermined problem, and therefore an 
infinite number of solutions exist. Under this circumstance, 
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the off-line algorithm must generate a solution (estimated 
health condition) “close enough” to the true solution; 
otherwise the algorithm may become unstable. In this paper, a 
simple constraint was imposed on the EKF: the estimated 
health parameters were constrained to change in the expected 
direction. The constraint was based on heuristic knowledge 
regarding health degradation, and it improved the stability of 
the EKF. Such an additional step is necessary for obtaining a 
reasonable solution in the underdetermined problem. 

Another way of improving the stability of the EKF is to 
utilize steady-state measurement data obtained at multiple 
operating conditions. At distinctly different operating 
conditions, the same health condition can induce different 
signatures in the measured engine outputs. The utilization of 
multi-point data augments the number of known parameters 
(sensor outputs) in the problem, and therefore can improve the 
stability and also estimation accuracy (refs. 19 and 20). 
Through this extra step, the robustness of the off-line 
algorithm may be improved, making the algorithm capable of 
running on an on-board engine computer. 

Conclusion 

A challenge in developing an on-line fault detection 
algorithm is making it adaptive to engine health degradation. 
If the algorithm has no adaptation capability, it will eventually 
lose its diagnostic effectiveness. To address this problem, the 
integration of on-line and off-line diagnostic algorithms was 
investigated. The off-line trend monitoring algorithm 
periodically estimates engine health condition over the course 
of engine’s life. Based on the estimated health condition, the 
on-line fault detection algorithm is updated. Through this 
integration, the on-line fault detection algorithm maintains its 
effectiveness while the engine continues to degrade over time. 

The integrated approach was investigated in a simulation 
environment using a nonlinear engine model. The evaluation 
result showed that this integration is essential to maintain on- 
line fault detection capability in the presence of health 
degradation. Practical issues of the integrated approach were 
also discussed. The communication between the on-line and 
off-line algorithm can be made less complex if the off-line 
algorithm is operated on an on-board engine computer. To 
make this possible, further research is needed in the area of 
off-line trend monitoring. 
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